Statistical fluctuations of pumping and rectification currents in quantum dots 
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We investigate the statistical fluctuations of currents in chaotic quantum dots induced by pumping 
and rectification at finite temperature and in the presence of dephasing. In open quantum dots, 
dc currents can be generated by the action of two equal-frequency ac gate voltages. The adiabatic 
regime occurs when the driving frequency is smaller than the electron inverse dwell time. Using 
numerical simulations complemented by semiclassical calculations, we consider both limits of small 
and large number of propagating channels in the leads when time-reversal symmetry is fully broken. 
We find that at intermediate temperature regimes, namely, ksT < A, where A is the mean single- 
particle level spacing, thermal smearing suppresses the current amplitude more effectively than 
dephasing. Motivated by recent theoretical and experimental works, we also study the statistics of 
rectified currents in the presence of a parallel, Zeeman splitting, magnetic field. 

PACS numbers: 73.21. La, 73.23.-b, 72.25.Dc 



I. INTRODUCTION 

The study of non-equilibrium transport phenomena in 
mesoscopic electronic systems has attracted much atten- 
tion in recent years, in particular for ballistic quantum 
dots formed in semiconductor heterostructures. Besides 
the investigation of non-linear I — V characteristics, it 
has been proposed that dc currents can be induced by 
the simultaneous application of two ac, shape-deforming 
gate voltages, the so-called pumping effectiMd An ex- 
perimental realization of an adiabatic quantum dot pump 
of electrons under these conditions has been reported; - 
For open quantum dots, where the system is connected 
to reservoirs by one or more propagating channels at 
each contact, the pumping current has usually both 
pure quantum (dissipationless) and rectified (dissipative) 
componentsi^ While attempts have been made to max- 
imize the former, it seems that current experiments are 
still dominated by the latter. Nevertheless, the effects of 
phase coherence can also be present in the rectified com- 
ponent. A clear indication of that has been given by the 
observation of a mesoscopic spin currenlii by Watson et 
al.^ in a setup where most likely rectification is still the 
primary source of pumping. 

In electronic mesoscopic systems, phase coherence and 
the consequent lack of self-averaging are responsible 
for the sample-to-sample fluctuations of most transport 
properties. For systems with an underlying chaotic dy- 
namics, these fluctuations are characterized by univer- 
sal statistics. Therefore, for either quantum pumping or 
rectification, it is important to understand what are the 
universal features of current fluctuations. Being univer- 
sal, these features can be described, in principle, by a 
theory based only on a few hypotheses about the sys- 
tem symmetries and its connection to the reservoirs and 



the environment. In fact, by comparing universal proba- 
bility distributions and correlation functions with actual 
experimental measurements, one may be able to gain in- 
formation about nonuniversal effects. Also, through this 
procedure, it might be possible to estimate the amplitude 
of phase decoherence present in the device. 

The theory of zero-temperature statistical properties of 
the pumping current in quantum dots^ explains qualita- 
tively several aspects of the experiment. - A quantitative 
approach, however, requires treating both temperature 
and the consequent dephasing effects. Previous work 
addressing dephasing did not account for the thermal 
rounding of the Fermi surface. ^'"'^^ Albeit the thermal con- 
volution for the pumping current being casted in a very 
simple expression, assessing the fluctuations requires the 
computation of nontrivial energy correlation functions. 
These correlation functions were calculated only in the 
absence of dephasing and when a large number of propa- 
gating channels {N ^ 1) is present . In this limiting 
case, the role of temperature is well understood. In the 
preset work we use analytical and numerical calculations 
to study the effects of temperature and dephasing in the 
case of small number of channels N > 2, which is more 
relevant to the experiments. Also, we revisit the case 
of iV 3> 1 from a semiclassical viewpoint, where results 
can be obtained in a more transparent and insightful way 
than from the diagrammatic technique. 

The paper is organized as follows. In Sec. ^ we briefly 
present the general theory of pumping currents in open 
QDs in the adiabatic regime. Finite-temperature pump- 
ing current fluctuations are expressed as integrals over 
correlation functions of scattering matrices in Sec. Ill Al 
In Sec. IIIBI we show that such correlation functions are 
readily calculated by means of the semiclassical approxi- 
mation. The small- A'^ limit is investigated numerically in 
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Sec. Ill CI The combined effect of dephasing and temper- 
ature is discussed in Sec. Ill Ul In Sec. Illll we address the 
influence of capacitive coupUngs on the pumping currents 
and their rectification effects. Mesoscopic fluctuations of 
spin pumping currents are discussed in Sec. IIVI Finally, 
in Sec. we present our conclusions. 



II. PUMPING CURRENTS 

We consider a quantum dot (QD) formed by later- 
ally confining a region of a two-dimensional electron 
gas (2DEG) through multiple electrostatic gate voltages. 
The electrons in the QD have free access to the rest of the 
2DEG through two point-contact leads. We assume that 
the "right" and "left" leads support Nji and Nl fully 
transmitting modes, respectively. The confining poten- 
tial of the QD undergoes a periodic shape deformation 
caused by the ac modulation of two gate voltages, which 
can be parameterized by Xi{t) = ai cos(wt + (f>i) and 
X2{t) = a2C0s{u}t + 02)- Here, cu denotes the driving 
frequency, ^ = 02 — </>i is the (constant) phase difference 
between voltages, and oi and a2 are their amplitudes. 
We focus our treatment of the problem to the adiabatic 
regime, when the driving frequency of the perturbations 
is much smaller than the electron escape rate (or inverse 
dwell time): u)<s:T/h = NA/2iTh, where N = N[i + Nl 
and A is the mean level spacing in the QD when isolated 
from the reservoirs. 

The adiabatic deformations of the QD confining poten- 
tial induce charge transfer through the point contacts. In 
the linear response limit, the amount of charge passing 
through a given contact / {I = L,R) is expressed as-^ 



(1) 



where the so-called emissivities can be written in terms 
of scattering matrix elements, namely. 



dni 
dX,, 



27r ^ ^ 



dX, 



(2) 



for i = 1,2. Thus, the total charge per cycle transfered 
through the contact I can be evaluated by integrating Eq. 
(|T|) over one period. 
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One can use Eqs. (O, (jSl, and rewrite the integral over 
time as an integral over the surface area A swept by Xi 
and X2 in parameter space over one period. Moreover, it 
is not difficult to prove that Ql = —Qr, since charge is 
not accumulated inside the QD (see below). As a result, 
we can cast the zero temperature pumping current going 
from left to right reservoirs as^ 



I,[Ep) ^"^^'£1 dX,dX2 K{EF,y^). 



(4) 



where X = (Xi, X2). Notice that only electrons in states 
at the Fermi energy contribute to the current (we have 
used the subscript in the current to indicate tempera- 
ture). In terms of the S matrix, the pumping response 



function IIq reads 



4(i?,X) = i^^Im 
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The energy and parameter dependences of the response 
function are determined by those of the scattering ma- 
trix. Therefore, it is necessary to make an explicit con- 
nection between the latter and the microscopic details of 
the system. Such connection is provided hyi^ 



1 



E-H{X) 



W, (6) 



where E is the energy of incoming and outgoing electrons 
(assuming that no bias voltage is applied to the system). 

is a Af X matrix that gives the coupling between the 
M resonant single-particle electronic states in the QD, 
described by the M x M Hamiltonian matrix -ff(X), and 
the N propagating modes in the leads. 

Throughout this paper we work with the hypothesis 
that the coupling between the QD and the leads is max- 
imal, meaning that there are no barriers for electrons 
entering and leaving the system. In this picture, W be- 
comes independent of any external parameter. This hy- 
pothesis is good inasmuch the changes in the confining 
potential due the periodic perturbation are local and take 
place far for the QD openings. This means that the dis- 
tance £ between the openings and the perturbed region 
must be such that £/Xf ^ 1. Should that not be the 
case, then dW/dX ^ 0, resulting in an additional pump- 
ing mechanism similar to that of a classical peristaltic 
pump, where transport occurs because the constriction 
are periodically opened and closed while the internal po- 
tential is varied. Although the existence of such classical 
effect in the experiments cannot be entirely ruled out, we 
do not take it into account in our analysis. The clear ex- 
istence of mesoscopic, sample-to-sample fluctuations and 
a marked temperature dependence of the magnitude of 
the current amplitude in the experiments seems to indi- 
cate that the peristaltic effect can be made rather weak. 
Thus, for a parametrically constant we have 



oXi oX, 



(7) 



where D ^ E - iJ(X) + iirWW'' . Using Eq. 0, it is 
straightforward to show that 11 = 11-^ = —11^. 

We are interested in chaotic quantum dots, where the 
Hamiltonian H can be modeled as a member of one of 
the Gaussian ensembles of random matricesii^ For this 
case, the Hamiltonian matrix elements are assumed un- 
correlated but equally distributed. Their variance A^/M 
determines the mean resonance spacing at the center of 
the band, A = nX/M . Several studies support that, in 
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general, the elements of dH/dX have themselves also 
Gaussian random entries (see, for instance, Ref. jlfl and 
references therein). We choose their variance to be 
(A/MXc)^iii In this way, we set the scale of X, making 
X = Xc correspond to the average parametric change 
necessary to cause one level crossing. In other words, 
Xi and X2 are measured in units of the average para- 
metric level crossing. Motivated by the fact that pump- 
ing currents are usually generated in the presence of an 
external perpendicular magnetic field that breaks time- 
reversal symmetry, in what follows we only address the 
case where H belongs to the unitary ensemble (GUE). 

The statistical theory presented above is an alterna- 
tive to the maximal entropy approach used to calculate 
the pumping current fluctuations in Rcfs. 3 10. At T = 
both approaches are equivalent. However, while the max- 
imum entropy has the advantage of leading to analytical 
expressions for the distribution of IIq at T = 0, it cannot 
be consistently extended to finite temperaturesiiS As we 
show bellow, temperature plays a very important role in 
suppressing the pumping current fluctuations. 



A. Thermal fluctuations 

Let us start addressing the thermal rounding of the 
Fermi surface, postponing the discussion of dephasing to 
Sec. Ill Dl The thermal smearing is easily accounted for 
by the integral 



/t(/x) 



dE 



dfr 

' dE 



Io{E), 



(8) 



where fr = {exp[(£^ — m)/^] + 1} (we assume fcs = 1 
hereafter). Hence, the pumping current variance reads 



(4> = 
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where Co is the pumping response autocorrelation func- 
tion at T = 0, defined by 



Co(£) = / dXi dX2 [ dYi dY2 D{e, X - Y) 

J A J A 



and 
i:>(e,X- Y) 



Ho(a* + |,x) Ho (/i-f,Y 



(10) 



(11) 



Here (• • •) indicates ensemble averaging. If the amplitude 
of the periodic perturbations is sufficiently weak, such 
that Xc, we can approximate the surface integrals 

by the mean value of the response function at X — 
times the area, yielding 



CQie)^A'^D{e,0). 



(12) 



At very low temperatures, the variance of the pumping 
current becomes proportional to Co(0), i.e., (Hq). In the 



absence of dephasing, we find that Hq has zero mean and 
variance given by (see Appendix 



var H 
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Note that var(Ho) diverges for N = 2. As N increases, 
this divergence is smeared out. For large N, Co(e) 
quickly tends to a universal, "semiclassical" form. The 
derivation of this function is presented below. 



B. Semiclassical pumping response autocorrelation 
function 



When the number of open channels N ^ 1, the 
pumping response autocorrelation function of Eq. H10|) 
can be directly evaluated using Miller's semiclassical S- 
matrixfiSi namely. 



Sa/siE, X) 



E 



p^(i;,X)e"^-(^'^)/'' , (14) 



where the classical trajectories that start at channel f3 
and end at channel a are labeled by ^{a <— P). Ac- 
cordingly, (T^, is the reduced action (with a Maslov phase 
included) and is the classical transition probability 
for going from /3 to a following the path /i. (Quantities 
indicated by a wide tilde are obtained in the semiclassi- 
cal approximation.) In the derivation of Eq. H14|l. the 
absence of tunneling barriers between the scattering and 
the asymptotic regions is implicit. 

Let us write the pumping response function as Ho = 
TT^^ X^aez The semiclassical approximation for 
Jq/3 reads 



j„^(i?,x) = J2 txidxl "^^ '''' 



(15) 



where fj, and v are trajectories connecting the same pair 
of two arbitrary channels a and (3. The parametric action 
derivatives are defined as^^ 



dX, 



, dH 
dt- — 

dX, 



(16) 



where the integral in evaluated along the trajectory /x 
over the time it spends in the QD. 

Our ergodic hypothesis is that averages over random 
matrix ensembles are equivalent to the energy aver- 
ages taken here. We average (/q^(_£/, X) over an en- 
ergy window 5E where the classical dynamics presents 
little changes, nonetheless fulfilling 5E ^ A. For this 
task, as customary, we neglect the energy dependence 
of the probabilities and keep only contributions from 
diagonal terms. This approximation is justified for tra- 
jectories with dwell times shorter than the Heisenberg 
time Tfj = h/A, since they are, in general, uncorrelated 
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for chaotic systems. Fortunately, without barriers, tra- 
jectories with t exceeding tjj are statistically negligible 
in the semiclassical regime oi N ^ 1. In the absence 
of system specific symmetries, the diagonal approxima- 
tion reads (exp[j(crp — o'i^)/Tt-]) se ~ ^a"'- implies that 
{Jc^piE, X))sE = and hence {Ilo{E, X))se - 0. 

We use a similar procedure to obtain the semiclassi- 
cal pumping response autocorrelation function D{e, 5X), 
defined as in Eq. Now we deal with a prod- 

uct of two Hq functions evaluated at different energies, 
E ± e/2, and parameter values, X ± 5X/2. The actions 
afj,{E±e/2, X±JX/2) are approximated in leading order 
of classical perturbation theory as 

a^{E±-,X±—)^a,{E,X)±t^-±Q,-—. (17) 

Let us start examining the J autocorrelation function. 
The diagonal approximation is used to compute the en- 
ergy average of sin[((Tp — a^)/h]sm[{afj_> — a^r)/h] and 
gives 



M,.'(q-/3) 



XyyPlJ.PlJ.'PiyPi''QlnQlfj.'Q2i^Q2i^' 

i 

h 



(18) 



X ex 



where J^^^ = JajsiE ± s/2,X± X/2). Since no special 
attention is payed to time-reversal symmetric paths, Eq. 
(|18|l represents the semiclassical correlation function for 
broken time-reversal symmetry. 

We proceed using classical sum rules to convert the 
sums in Eq. (|18|l into time integrals. For that purpose, or- 
bits are grouped with respect to common traversal times 
and averages are taken within these sets. In what fol- 
lows, we describe the details how this procedure is im- 
plemented following two basic steps. 

For chaotic systems the transition probabilities fol- 
low the analogue of the Hannay-Ozorio de Almeida sum 
rule for open systemsf^ 



J2 

t<tu<t+5t 



= ^e~'/^St^p{t)St, (19) 



where J2t<t, 



<t+st 



is the sum of all classical transition 



probabilities following the trajectories belonging to a 
time interval [t,t -f 6t], where 5t is classically small. In 
Eq. l(T^ . the decay time is r = h/T^, where F is the 
escape width, also known in this context as the Weisskopf 
S-matrix autocorrelation length. 

The transition probabilities and the parametric ac- 
tion derivatives are uncorrelated. In addition, Qi^ 
and (52/j are uncorrelated, provided /i dwells inside the 
QD for a couple traversal times. Thus, for fixed X, 
the time-average in Eq. H19|) runs over a large number 



of scattering orbits with Qp(t), which may be consid- 
ered as samples of the probability distribution Pf(Q) = 
Pt{Qi)Pt{Q2)- The latter is assumed to be Gaussian, 



Pt{Q^) 



1 



: exp 



2Qf(i) 



(20) 



The Gaussian width is a function of X and E, known to 
grow diffusively with time,^^ namely. 



1 



(21) 



t<t^<t+st 



where Af{t) is the number of trajectories /i within the 
time window [t,t + dt\. The diffusion constant B is not 
expected to depend on i since both Xi and X2 are generic 
shape parameter. (This is not the case if, for instance, 
Xi would stand for a magnetic fiend and X2 for a shape 
deformation.) B is given by 



B{E, X)^2j^dt — (p(t), q(i)) — (p(0), q(0) 

(22) 

It is important to notice that the decay of the classical 
correlation function in Eq. H22|l need only be integrable 
(many available chaotic systems do not exhibit full expo- 
nential decay of the correlations). A detailed discussion 
of the classical properties of and its semiclassical im- 
plications for density correlation function can be found 
in Ref. 20. In particular, there it is shown that Xc is inti- 
mately related to the classical diffusion constant, namely, 

X,^\h^/Bf/yTT^ 

By taking the Gaussian average over Qi^, inserting 
(|19|l and H22I) into (|18|) . and integrating over time, we 
arrive at 



D{e,6X) 



var n 







with 



varllo 



2(<5X)^ 



1 



(23) 
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(24) 



This expression agrees with the random matrix theory re- 
sults of Ref. J^, also obtained for ^ 1. The semiclas- 
sical var Ho is also consistent with Eq. (|13l) . as it should. 



C. Small- Af limit 

It is a difficult technical task to calculate an analyt- 
ical expression for Co(£) in the small- A^ limit, which is 
the case of experimental interest. To bridge this gap, 
we relied on numerical simulations. We typically gener- 
ated 10^ members of an ensemble of S matrices defined 
by Eq. ©, following the prescription given in Ref. m 
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FIG. 1: Pumping current autocorrelation function 
Co{e)/Co{Q) versus energy e in units of F. The solid 
circles represent the A'' = Nl + A''^: = 4 case, whereas empty 
circles, triangle, and crosses stand for A'' = 6, 8, and 10, 
respectively. In all cases Nl ~ Nr. The solid line is the 
autocorrelation function given by Eq. II23II . 



FIG. 2: Pumping current variance {l'^) in units of (eti;^/27r)^ 
versus temperature T in units of A. The thick solid line stands 
for the A'^ = Nl + Nr = 4 case, whereas dashed, dotted, and 
short-dashed lines stand for A'^ = 6, 8, and 10, respectively. In 
all cases Nl = Nr. Inset: Same cases, but for the pumping 
current variance scaled by Co(0) (temperature measured in 
units of F). 



The results for > 2 are illustrated by Fig. ^ No- 
tice that as N increases, the semiclassical limit is at- 
tained very fast. Nevertheless, for N = 4, which is the 
experimental case in Ref. |^ the depart from the semi- 
classical limit is still significant. The effect of having a 
small iV is best captured in Fig. |21 where we show the 
pumping current variance in units of {eujA/2Tr)'^, namely, 
(Z^) = {It'^) /{eujA/2TTy , as a function of temperature. 
The numerical results were obtained by carrying out the 
integration in Eq. through an interpolating of the 
data shown in Fig. ^ We can observe that the fluctu- 
ations decrease with increasing temperature and N. In 
the inset we scaled (I^) by C(0) and the temperature by 
F. We observe that while the curves for iV > 8 tends 
to fall onto each other, the = 4 case shows a marked 
different behavior. 

The limit case of iV = 2 is special: the variance of the 
pumping response function Uq diverges due to the long 
tails in the distribution IIo;- Hence, in the absence of 
dephasing processes, (/|,) also diverges. We address this 
issue in the following subsection. 



D. Dephasing 

The origin of the quantum pumping current is inter- 
ference. Thus, the effect is susceptible to dephasing cre- 
ated by the interaction of electrons with phonons, pho- 
tons, and fluctuations in the electromagnetic environ- 
ment. While a precise, accessible microscopic theory of 
dephasing for open quantum pumps is still lacking, some 
quantitative results can be obtained through phenomeno- 



logical modelsAiS In particular, the voltage probe model 
of Biittiker— provides a simple way of introducing de- 
phasing by adding a third lead to the QD which neither 
inject nor drain a net current. Electrons can move in and 
out of the QD through this lead and loose phase coher- 
ence in the process. The amount of dephasing can then 
be tuned by changing the characteristics of the third lead, 
such as coupling constant p and number of channels A^^. 
It is customary to use a single parameter = pN^ to 
parametrize the dephasing. We take ^ 1 and p <^ 1, 
while keeping constant. 

Here we adopt the formulation of Ref. and break up 
the contributions to the pumping current into two parts, 



(25) 



The separation into two components is convenient be- 
cause Hq"^ becomes Hq [Eq. 0] as the coupling between 
the third lead and the QD goes to zero (see below) , while 
Hq'^'^ disappears in the same limit. The latter represents 
the contribution to the pumping current coming from a 
voltage rectification effect (See Sec. IIII|) . 

Using that Hq + H^ + Hq = to enforce current con- 
servation (where </> denotes the third lead), one finds the 
relationiS 



h; 



G 



R4> 



H 



Gr 



dir 



(26) 



where G^^ and Gfj^ are the conductances of the QD 
between the third lead and the left or right leads, respec- 
tively. Through steps analogous to those taken to write 
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FIG. 3: Pumping response variance C^iiQ) = var(n^) as a 
function of P^. Solid circles are the result of our numerical 
simulation (solid lines serve as a guide to the eye), (a) A'^ = 2 
case. Inset: Same as before, but magnifying the small 
range (statistical standard deviations shown), (b) TV = 4 
case. 



Eqs. to one finds that the rectification current is 
given by 



\ '^^ dXi J dX2 Gl4, + Gr^ 

(27) 

Below, we study the statistical fluctuations of IIw 
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In Fig. 131 we show the effect of dephasing in the pump- 
ing response variance, C^{0) = var(n0), as a function of 
the dephasing parameter for = 2 and 4. When 
N = 2, the long tails in the distribution of Hq make the 
numerical assessment of var(n0) increasingly difficult as 
goes to zero (see inset of Fig. ■ To investigate this 
problem, we generated 10^ realizations of the S matrix 
for a range of energies E and P^ < 0.1. We found that the 
standard deviation increases with the number of realiza- 
tions; moreover, large values of jll^l were accompanied 
by abrupt fluctuations of II^ with E. Upon shrinking 
the energy steps in our numerical calculations, we found 
that the number of large fluctuations decreased, but their 
amplitudes increased. This scaling procedure was com- 
putationally very costly; yet, despite the somewhat poor 
statistics, our results suggest a fractal behavior of II^ as 
a function of energy, typical of power-law distributionsi^^ 
Therefore, our simulations strongly support that C^{0) 
diverges for = 2 as P^ — > 0, in agreement with Ref. 
U The apparent fractal nature of I1^{E) in this case in- 
dicates that exchanging the order between thermal and 
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FIG. 4: Normalized pumping response II^ autocorrelation 
function C^(e)/C^(0) versus e/F^ for Nr = Nl = 2. The 
solid circles, squares, triangles, empty circles, and crosses 
stand for P4, = 0.25, 0.5, 1.0, 2.0, and 3.0, respectively. The 
solid line stands for Eq. I23II . 



ensemble averages would likely not render a finite value 
for {I^) either. 

For larger values of the dephasing parameter, say 
P<P ^ 0.1, we obtained a converged value var(n0) (small 
standard deviation) with 10^ realizations, despite of the 
very large values of the higher moments. For P^ > 1, 
convergence was attained already with 10* realizations. 
We observed (not shown here) that our numerical simula- 
tions a gree with the analytical results for C^{0) presented 
in Ref.lin. 

In practice, at very low temperatures, var(n0) is ex- 
tremely sensitive to small changes in P^. In this case, 
although in principle possible, it is very difficult to give a 
quantitative numerical description of both dephasing and 
temperature effects. We do not pursue this path here and 
direct the discussion instead to > 2, where certainly 
no divergences occur. 

A very strong dephasing dependence is absent for 
N = A. Here, our simulations show that C^{Q) remains 
finite as P^ approaches zero, as depicted by Fig.|3jD and in 
quantitative agreement with Eq. H13|) . In Fig. ^ we show 
the effect dephasing in the dimensionless pumping energy 
autocorrelation function. The correlation function (e) 
is defined as in Eqs. (HU and We Nr = Nl = 2 
and vary the dephasing parameter P^ . The results repre- 
sent an average over lO'' realizations. We find empirically 
that the energy correlation length scales as 



(28) 



Similarly to the case without dephasing, as P^ increases 
one quickly reaches the semiclassical regime, character- 
ized by the universal correlation function given by Eq. 
(|23|l . We verified (not shown here) that our simulations 
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1.0 1.0 



FIG. 5: (X^) for TV = 4 as a function of the temperature and 
the dephasing strength. 



for Cd,(0) coincide with the analytical results obtained in 
Ref. 10, as expected. 

At this point we can compare the relative role of 
dephasing and temperature for the pumping current. 
This is done in Fig. O where we display (I^) = 
{lT'^)/{eujA/2TTy as a function of T and P^. For Nr = 
Nl > 2 we find that, at low temperatures, where typi- 
cal QDs have for instance T/A w 0.5 and « 0.5^^ 
temperature plays a significantly more important role 
than dephasing. The exception is the remarkable case 
of Nfi = Nl = 1 where, in the absence of dephasing, the 
variance of the pumping current diverges. 



conductances to vary in time, an incoherent dc current 
can flow between left and right reservoirs, as long as Gl^ 
and Gr^ are nonzero. This effect is intrinsic to the QD 
and has therefore been included in the expression for the 
quantum pumping current. We remark that the rectified 
current in this case is also asymmetric with respect to 
magnetic field inversion.^ft For an open QD, this contri- 
bution to the current vanishes as the dephasing strength 
decreases. 

Secondly, rectification also appears due to the capac- 
itive coupling between the shape-deforming electrodes 
and the left and right reservoirs.^-'^- The displacement 
currents and the conductance of the QD oscillate with 
the same frequency, producing a net dc charge current 
between the reservoirs. The exact way by which the dc 
current is induced depends on the particular measure- 
ment setup. To illustrate this point, in Fig. |Blwe show 
schematically the two equivalent circuits, namely, for (a) 
voltage or (b) current measurements.^* Let us briefly de- 
rive expressions for the QD voltage and current for these 
setups, expanding the discussion found in Ref. |^ 
(a) Voltage setup: in this case the reservoir-QD-reservoir 
loop is open. Calling V = Vr — Vl, with Vl = 0, we find 
that 

I^C2R-^^{X2-V) + CiR-^^{X,-V). (29) 
Since / = G V" as well, we arrive at 



V = 



G.4(x.-y) + G,4(x,-y) 



(30) 



For sufficiently adiabatic pumping (MHz) and small ca- 
pacitive couplings (pF), we can have V <^ ^1.2, given 
that G-^ <26Kn. As a result,^ 



III. RECTIFIED CURRENTS 

In AlGaAs/GaAs quantum dots formed by lateral elec- 
trostatic gates, the rectified component of the dc cur- 
rent tends to dominate over the quantum pumping one 
at low driving frequencies. This has been recently veri- 
fied in an experiment where the symmetry properties of 
induced currents with respect to an external perpendicu- 
lar magnetic field were studied. It was found that while 
pure quantum pumping currents should be asymmetric 
with respect to field inversion, the actual measured 
current showed a strong even symmetry, characteristic 
of rectification.^ It is therefore important to characterize 
the statistical fluctuations of the rectified current as well. 

Rectification can arise from two different sources. 
Firstly, as mentioned in Sec. IIIDI it appears when in- 
elastic processes take place inside the QD, leading to de- 
phasing. When one models the dephasing process by 
allowing carriers to move in and out of the QD incoher- 
ently through a third lead, an additional reservoir with 
time-dependent chemical potential, /^^(t), is required^ 
Since the time modulation of the QD shape causes the 



V 



1 fr^ dX2 

G r^^^ 



Gi. 



dXi 

'~dr 



(31) 



(b) Current setup: in this case both reservoirs are 
grounded and Vr = 0. Since ig = I — iiL — i2L, we 
find 

I=^ + C2Lj^{V-X2) + C^Lj^{V^X,). (32) 



Using I — GV, we get 



RG 
l-RG 



C2LJ^{X^ 



V) 



Gii-(X2 



V) 



(33) 

Here we can also neglect V with respect to Xi^2 on the 
right-hand side of Eq. (|33|) . Moreover, in the MHz 
range, where the pumping experiments are carried out, 
the impedance of the current meter should be smaller 
than the quantum dot resistance, namely, RG ^ 1^ 
Thus, we finally obtain^ 



I^RG\C2l'^-^ + cJ^' 



dt 



dt 



(34) 
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FIG. 6: Equivalent circuits for the quantum pumping mea- 
surements: (a) voltage setup; (b) current setup. The elec- 
trodes are coupled to the right (left) reservoir through ca- 
pacitances CiR and C2R {C'lL and C2l), with displacement 
currents im and i2R {iiL and 22^). In (a) the reservoir- 
dot-reservoir circuit is open, while in (b) both reservoirs are 
grounded. R accounts for the circuit and the current meter 
internal resistance. Xi and X2 are the ac voltages in the 
electrodes and G{t) is the quantum dot conductance. 

Focusing now in the current setup [Eq. H34|) ]. it is 
straightforward to show that an expression for the recti- 
fied current similar to Eq. Q exists, with the kernel 



n[;"p = Rc, 



dG 



RCi 



dG 
dx'2 



(35) 



where we call dG/dXi the parametric conductance ve- 
locity with respect to Xi. The statistical fluctuations 
of this quantity were investigated in Ref. l30l where the 
distributions of dG/dX were presented for the N = 2 
case. Recall that the linear conductance G is given by 
the Landauer formula 

c2 



G 



90, 



(36) 



with the dimensionless conductance at zero temperature 
defined as 



90 



E 

iieR 



/3l 



(37) 




1.0 2.0 

e /r 



3.0 



FIG. 7; Parametric conductance velocities autocorrelation 
function Cggg/gxis) as a function ofe/T. The solid line stands 
for Eq. II42|I . while squares, circles, and triangles stand for 
Nl ~ Nr — 1, 2, and 3 respectively. 



It is difficult to estimate the amplitude of the rectified 
current derived from XIq^^^, since it depends on R which is 
unknown. However, the symmetries and statistical prop- 
erties can be determined and should differ from those 
of the pure quantum pumping current component. For 
instance, notice that XIq^^ depends directly on dgo/dX. 
Since the conductance is symmetric with respect to field 
inversion, and so will be the rectified component of the 
current in Eq. H34|) .^ Let us focus now on the statisti- 
cal properties of IIq^'' and on its temperature and de- 
phasing dependences. We argue that thermal smearing 
plays a key role in suppressing fiuctuations of the recti- 
fied current as well. In App. |21 we study a very closely 
related issue and show how thermal smearing explains 
the discrepancy between the recently experimental mea- 
sured parametric conductance velocity distribution^ and 
the analytical predictions at zero temperature^ 

A simple inspection of Eqs. © and {Tj) shows that 
(dgo/dXi) ~ 0, since it is proportional to {dH/dX. 
Hence, (Hq^'P) =0 
hand, is given by 



The variance of Hq^^, on the other 



/ Re^ 

var(Hr) = ( ^ 





^^'^A-dx.-dx, 



(38) 



For chaotic systems, provided the two parametric per- 
turbations X\(i}j and Xi(t) have the same amplitude, we 
expect that 



(39) 




9 



If the periodic perturbations Xi and X2 are acting 
at locations far (by several Fermi wavelengths) from 
each other, it justifiable to assume that the paramet- 
ric conductance velocities with respect to different Xi 
are uncorrelated, {{dgo/ dXi){dgo/ 8X2)) = 0, leading to 
var(n;5''P) = R^{Cl+Cl){{dgQ/dXf). This assumption 
leads to an analytical expression for the zero-temperature 
rectified current variance, var(i''°'^*), in analogy with Eq. 
Q, since 




\NlNl 



(40) 



[Details on the derivation of Eq. 1)40(1 are presented in Ap- 
pendix^] In order to compute the temperature depen- 
dence of var(i'''"^*) we need first to calculate the paramet- 
ric conductance velocity energy autocorrelation function, 
namely. 



2J axr"2 



(41) 



(recall that {dgo/dX) —Q). In Fig. [7|we show the results 
of our simulations for various values of N . Again, as N 
increases Cgg^/gxis) very rapidly converges to 



Cdg„/dxi£) - 1 +'(e/r)2 ' ^^^^ 

where Cggg/gxi^) is given by Eq. (^UJ), while deviations 
are quite large for small N. As in Sec. Ill CI the thermal 
fluctuations are enhanced in the small-channel case. 

For completeness, let us discuss now the oppo- 
site, correlated case, where {{dgo/dXi){dgo/dX2)) — 
[dg^/dXf. Nowvar(n^^P) = R^{Cr-C2?{{dgo/dXf). 
making and the effect of rectification very small when 
the capacitive coupling is close to symmetric (Ci ~ C2). 
However, it is unlikely that such condition is satisfied in 
real experimental setups. 

Dephasing effects can also be included in IIq'^'' phe- 
nomenologically by using the voltage probe model in its 
original formi^ For that purpose, the dimensionless con- 
ductance is replaced by 

5 = 90 + • (43) 

9l4> + gB4 

We have calculated the distribution of dg^/dX for sev- 
eral values of N and the dephasing parameter P^, as 
introduced in Sec. IIIDI The conductance velocity auto- 
correlation function scales in the same way as in the case 
without dephasing, provided we replace F by F^. Results 
for the case N — 2 are shown in Fig. |S1 The dependence 
with is weaker than in the pure quantum pumping 
case (compare with Fig. 0}. 



IV. SPIN CURRENTS DUE TO PUMPING AND 
RECTIFICATION 

It was recently pointed out that it is possible to pump 
a spin current without any net flow of charge through a 




FIG. 8; Parametric conductance velocities autocorrelation 
function Cgg^/gxis) as a function of e/T^ in the presence 
of dephasing for Nl = Nr = 1. The solid line stands for Eq. 
14211 : squares, circles, triangles, and inverted triangles stand 
for =0.5, 1, 2, and 4, respectively. 



QDii The basic idea is to use a parallel magnetic field 
to Zeeman split resonant states around the Fermi en- 
ergy. By doing so and tuning the QD shape, one can find 
a situation where up and down spin components of the 
pumping current flow in opposite directions. The mech- 
anism by which charge is pumped can be either the ideal 
quantum pumping of Sec. ^ or the rectification due to 
the capacitive coupling between pumping electrodes and 
reservoirs, as described in Sec. IIIII A recent experiment 
has observed the effect 1^ 

An important practical question is by how much is 
the spin current attenuated by dephasing in the orbital 
(charge) sector of the wavefunctions. Furthermore, it is 
also interesting to know the amplitude of the effect when 
the spin pumping current is caused by rectification alone 
and the number of propagating channels in the leads is 
not large (both cases had not being considered in the 
analysis of Ref. (3) 

The dependence of the spin current polarization am- 
plitude on the applied parallel field and temperature can 
be estimated from the correlator (/|/|). Here /|.| = 
/T(Ai ± Ez/2), where Ez = g*^B-B||/2 is the Zeeman 
energy. Using the relations presented in Sec. m it is 
straightforward to show that the correlator can be writ- 
ten as 

/■°° (7 r / ^ XI -2 

{I^I^) = T'J^ de— [2Tsinh(— jj Coie + Ez). 

(44) 

This expression is valid regardless of the nature of the 
pumping mechanism. 

As discussed in Sec. IIIBI in the case of pure quantum 
pumping, Eq. (|23|l provides a semiclassical approxima- 
tion to the energy correlator in the absence of dephasing 
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FIG. 9; Spin pumping current polarization coefficient as a 
function of Zeeman energy Ez (parallel magnetic field), (a) 
Curves for different temperatures and a fixed number of chan- 
nels {Nr = Nl = 2). (b) Dependence on the number of 
channels for a fixed temperature. 



and when N ^ 1. Using that expression and carrying 
out simple manipulations, we find that Eq. H44f) becomes 
identical to Eq. (3) of Ref. 0, namely, 



/•OO 

(/t/i) = Co(o)r / dT(i + rT)e-r^ 

Jo 



ttTt 



sinh(7rTT) 



(45) 



The latter was originally derived in Ref. from the 
linear-response limit of the large- iV diagrammatic formal- 
ism developed in Ref. Notice that when dcphasing is 
present, F has to be replace by F^, as defined in Eq. (|2Hll, 
but the functional dependence on T and Ez remains the 
same. 

In the small- TV limit, we do not have a closed analytical 
expression for the energy correlator. In this case, the 
dependence of (Iili) on temperature and magnetic field 
amplitude can only be obtained numerically. Our results 
for this case are presented in Fig. |51 together with curves 
derived from Eq. H45I) . We have opted for plotting the 
spin current amplitude in terms of the spin polarization 
coefficient, defined as 



1 



(46) 



Notice the rather weak temperature dependence after 
rescaling the Zeeman energy. In contrast, the dependence 
on the number of channels is much more pronounced; 
consequently, a strong dependence on the dephasing pa- 
rameter Pff, also occurs. Thus, the same restrictions to 
orbital (charge) quantum pumping due to decoherence 
by the environment apply to the spin case. 
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FIG. 10: Same as in Fig. El but for rectification-induced spin 
pumping currents. 



Equation (|44|l is also valid when rectification is present 
and dominates the pumping mechanism. In this case, in 
the large- limit, the energy correlator of Eq. (|42|l can 
be used, leading to an expression very similar to Eq. 145|l . 
namely, 

cos{E;,t). 
(47) 

As in the dissipationless case, here we have to rely on nu- 
merical calculations to obtain the correlator in the small- 
N limit (see Sec. Illl|) . Comparative results are shown in 
Fig. HlOfl as a function of Zeeman energy (magnetic field) 
and for different temperatures and number of channels. 
The curves are very similar to those of Fig. |^ The depen- 
dence on temperature is just slightly more pronounced for 
capacitive pumping. Considering that the energy corre- 
lators for the two pumping mechanisms are very similar 
(a Lorcntzian square for pure quantum as opposed to a 
simple Lorentzian for the capacitive case, in the semiclas- 
sical limit), this result comes as now surprise. We con- 
clude that temperature and dephasing dependences alone 
are not strong indicators of the nature of the pumping 
mechanism, both for the charge and the spin cases. 



V. CONCLUSIONS 

In this work we presented numerical (small-A^) and 
semiclassical (large-iV) calculations of statistical mea- 
sures of charge and spin pumping currents in open chaotic 
quantum dots. Both pure quantum and rectification 
pumping mechanisms were considered when time-reversal 
symmetry is broken and the regime is adiabatic. We paid 
particular attention to the dependences on temperature 
T, number of propagating channels N , and dephasing. 
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We were able to draw several conclusions from our re- 
sults. 

Our initial motivation was to propose measuring the 
damping of pumping current mesoscopic fluctuations as 
a way to estimate the amount of dephasing present in 
the experimental setups. However, we found that, for 
realistic conditions, thermal smearing tends to be a more 
important effect than dephasing. Lower temperatures are 
required (T <^ A) for dephasing to become quantitative 
and qualitative distinguishable. 

Perhaps one situation where dephasing does have a 
stronger impact than thermal smearing is when N = 2, 
namely, when there is one single propagating channel in 
each lead. In this case, dephasing cuts off the divergence 
of the zero-temperature pure pumping current variance. 
Taking the thermal average before the ensemble average, 
on the other hand, does not seem to yield a finite vari- 
ance when dephasing is absent. In fact, our numerical 
simulations suggest that the zero-temperature pumping 
current response function, Hq, has a fractal behavior for 
N = 2. A long tail in its probability distribution was 
observed. Events belonging to this tail were connected 
to high but isolated peaks in Hq as the energy is varied. 
We observed that such peaks proliferate as the energy 
increment decreases. Albeit not yet accessible to experi- 
mental investigation, a more quantitative understanding 
of the suppression of this fractal pattern by dephasing is 
necessary. 

We also derived an expression for the pumping re- 
sponse autocorrelation function using the semiclassical 
approximation in the large- A'^ limit. Our result coincides 
in the leading order in powers of 1 /N with that reported 
m Ref. in which was obtained by the diagrammatic 
technique. The simplicity of the semiclassical derivation 
relating the pumping current to the instability of the un- 
derlying classical orbits provides additional insight to the 
physical process. In addition, our simulations show that 
the semiclassical limit is very quickly attained. We also 
calculated the pumping response variance using the S- 
matrix maximum entropy approach for pure quantum 
and rectification pumping. The results are valid for any 
number of channels and are both consistent with the nu- 
merical simulations and with the results presented in Ref. 

M 

For the rectification case, we found that the expres- 
sions relating the pumping currents and voltages to the 
driving perturbations are somewhat more involved than 
that proposed in Ref. 0. Our results only coincide with 
those of Ref. in the limit of weak capacitive coupling 
between plunger gates and leads, and when frequencies 
are not too high (so that the internal impedance of the 
current meter remains lower than the quantum dot resis- 
tance). 

Rectification currents induced by capacitive couplding 
also show mesoscopic fluctuations. However, dephasing 
effects are less pronounced than for pure quantum pump- 
ing. Moreover, the variance does not diverge at zero tem- 
peratures for any value of A^ > 2. 



Pure quantum pumping spin currents have the same 
characteristics of the charge case. Namely, their vari- 
ance diverges in the absence of dephasing. We found 
that increases in temperature, dephasing, or A'^ suppress 
spin polarization in a similar (strong) way. New mecha- 
nisms of pure spin currents generations based on pump- 
ing have been proposed recently. Particularly attractive 
are those based on the spin-orbit coupling present in two- 
dimensional electron gases formed in III-V heterostruc- 
tures, which do not require the application of large mag- 
netic fieldsj^^*^ For the future, we plan to study how 
thermal smearing and orbital and spin dephasing affect 
the magnitude of these new mechanisms. 
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APPENDIX A: MEAN AND VARIANCE OF THE 
PUMPING RESPONSE FUNCTION 

Here we derive Eq. H13|) . It is convenient to 
parametrize the S matrix, its parametric-derivatives, and 
its energy derivative as 



s = uv, 



dS 
dx' 



B S 2ni 
lUQx.V, and — = —UQeV, 



dE A 

(Al) 

respectively, for j — 1,2. Here, U and V are N x N 
unitary matrices uniformly distributed over the unitary 
group and independent of Qxj and Qe- In turn, the 
latter are N x N Hcrmitian matrices satisfying the joint 
distributioniS 



P{S, Qe,Qx„Qx,) cx (det Qe) 



-9N/2 



X exp 



-tr 



1 ^ 



Also, we parametrize Qxj as 

Qx, - ^^''k,^-\ 
where ^ is a complex N x N matrix such that 



Qi 



(A2) 



(A3) 



(A4) 



and Kj is an Hermitian N x N matrix whose elements 
are Gaussian distributed with zero mean and variance, 

{KabKa'h') = 4:{Sab'6ba' + Saa'Sbb'), (A5) 
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FIG. 11: Distribution of parametric conductance velocities 
P{dgo/dX) for N = 2. The histogram is our numerical sim- 
ulation, whereas the solid line stands for the analytical result 
of Ref. fSdL In the inset we show the tails of the distribution 
in a log-log scale. 



as can be easily seen substituting Eqs. ljA3|l and (|A4|) 
into the distribution (|A2I) . 

We substitute Eq. ljAl|l into Eq. |SJ| and average over 
the matrices [/, V using the results of Ref. |^ Next, we 
average over Qx^ using the parametrization of Eq. (|A3() 
and Eq. (|A5p . As a result, we obtain (Ho) = and 



FIG. 12: Distribution of parametric conductance velocities 
Pidgr/dX), Nr = Nl = 1, and T = 0.6A. The solid line 
indicates the best Gaussian fit. 



APPENDIX B: PARAMETRIC TRANSMISSION 
COEFFICIENT VELOCITY VARIANCE 

In this appendix we show how to obtain {{dga/dX)'^), 
as given by Eq. (|40l) . As in Appendix 1X1 we find useful 
to parametrize the S matrix as in Eq. (|A1() . In this case, 
Qx, Qe, and S have the joint distribution^Si 



var Ho 



SNlNr 
' N{N^ - 1) 

N 



N 



a.b.c—1 



((Q|)aa(Q_E)b6(Q£;)cc) 
(A6) 



Now, we write in its diagonal form as Qe = AtA^ , 
where f is the eigenvalue matrix and A is a random uni- 
tary matrix. Substitution of the diagonal form of Qe 
into Eq. l|A6p gives a result independent of A such that 
the average over A is easily done. After some algebraic 
simplifications, we arrive at 



varfn 



0) 



32NlNr 

N + 1 
+ {N-2) (riVarg) 



(rM)+2(rfr,) 



(A7) 



Finally, since the dimcnsionlcss scape rates a;„ — 1/t„ 
are distributed according to the generalized Laguerre 
ensemblei^ 



P(xi, 



7i<m 



X„ 6 



(A8) 



The averages appearing in Eq. IjAVp can be calculated 
by direct integration. The result is Eq. H13|) . 



P{S,Qe,Qx) oc (det Qe 



-4JV 



X exp <— tr 



Qi 



8 



{Qe'Q^ 



(Bl) 



We use Eq. IjAip to express {{dgo/dX)^) as a function of 
U, V, Qx, and Qe- With the help of Ref. l3ll we average 
over the matrices U and V. Next, using Eqs. ljA3|) and 
^5|). we integrate over Qx to obtain 



dgo 
dX 



Nt,N, 



N{N^ - 1) 

N 



Re 



N 



a, 6=1 



bb/ 



(B2) 



Again, we write Qe in its diagonal form as Qe = AtA^ 
and substitute into Eq. (|B2|I . After some algebra, we 
arrive at 

\[dx) / = NiN-l)iN+ir 

. . 

The averages over and tiT2 can be done explicitly by 
direct integration using Eq. IjASp . The result is Eq. H4U|I . 
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APPENDIX C: PARAMETRIC CONDUCTANCE 
VELOCITY DISTRIBUTIONS 

This appendix serves to remedy the discrepancy be- 
tween the theoretical^S and experimental— distributions 
of parametric dimensionless velocities dgo/dX. This dis- 
cussion gives further support to the statistical theory em- 
ployed in this paper, and provides a further illustration of 
the important role of temperature in the statistical fluc- 
tuations of conductance-related quantities. In addition, 
we illustrate the accuracy of the simulations by compar- 
ing our results to the theoretical P{dgo/dX). 

For N = 2 the distribution of dgo / dX is known for all 
symmetry classes It shows a singularity at zero deriva- 
tive: a logarithmic divergence in the presence of time- 
reversal symmetry (/3 = 1) and a cusp in the absence 
of that symmetry {f3 — 2). The tails of P{dgo/dX) are 
algebraic and follow 

P{dgo/dX) cx {dgo/dX)-^-^. (CI) 



In Fig. II II we show a comparison between our numerical 
results and the analytical distributions of dgo/dX. 



By numerically convoluting go{e) with the thermal dis- 
tribution, we obtain gT- The histogram is show in Fig. 
^] for T/A — 0.6, corresponding to the experimental 
situation in Ref. |2y. Notice that even at such low tem- 
peratures, T is slightly larger than the dgo/dX energy 
autocorrelation length, suppressing large fluctuations of 
go at zero temperature and favoring a Gaussian distribu- 
tion, as observed in the experiment. 



The quality of this result suggests that a direct com- 
parison between experimental and theoretical values of 
Yai{dgT I dX) can be used to set the scale of the X pa- 
rameter. The distribution of rectified currents 
can be obtained straightforwardly from P{dG / dX). 
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